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Abstract 

We propose a scenario of quantum memory for light based on Raman scattering. The storage 
medium is a vapor and the different spectral components of the incoming signal are stored in 
different atomic velocity classes. One uses appropriate pulses to reverse the resulting Doppler phase 
shift and to regenerate the signal, without distortion, in the backward direction. The different 
stages of the protocol are detailed and the recovery efficiency is calculated in the semi-classical 
picture. Since the memory bandwidth is determined by the Raman transition Doppler width, 
it can be adjusted by changing the angle of the signal and control beams. The optical depth 
also depends on the beam angle. As a consequence the available optical depth can be optimized, 
depending on the needed bandwidth. The predicted recovery efficiency is close to 100% for large 
optical depth. 

PACS numbers: 42.50.Ex,42.50.Md,03.67.-a 
Keywords: quantum information, storage 
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I. INTRODUCTION 



The storage of quantum information in an atomic ensemble has received a great deal 
of attention over the past ten years or so. Protocols based on electromagnetic induced 
transparency (EIT) have been investigated by many groups both theoretically 1, J and 
experimentally, leading to storage and retrieval demonstration of both discrete % |5( and 
continuous js, 7, ^] quantum variables. Although successful, the EIT storage scheme suf- 
fers from time-bandwidth product limitations. In EIT-based protocols, the reduction of 
group velocity is used to spatially confine the input signal within the boundaries of the 
storage medium. Simultaneously, the signal spectrum must not exceed the bandwidth of 
the transparency window associated with EIT. Since the group velocity is inversely pro- 
portional to the width of the transparency window, it follows that, the larger the storage 
bandwidth, the shorter the temporal profile the memory can accommodate. This protocol 
has been demonstrated using level schemes in which inhomogeneous broadening does not 
play a significant role. Extension to inhomogeneously broadened systems does not improve 
the time-bandwidth product capabilities. 

On the other hand, inhomogeneous broadening can be of critical importance in other 
protocols for storing quantum information. For example, Doppler broadening determines 
;he storage bandwidth when a signal pulse is totally absorbed in an optically dense medium 
9|. While the temporal components of the input signal are distributed in atomic state 
coherence along the axial direction in EIT, the spectral components of the input signal 
are spread over the inhomogeneous frequency distribution of the atoms in the absorption 
protocol. Just as in EIT, information is stored in a long lifetime Raman coherence in the 
absorption protocol, but the maximum duration of the input signal is independent of the 
storage bandwidth, being ultimately limited only by the inverse homogeneous line width. 
The resulting time-bandwidth product capacity, given by the ratio of the inhomogeneous and 
homogeneous widths, is reminiscent of the photon-echo based storage techniques that were 
developed in the past 10J. However, unlike some of those classical light storage schemes, 
the proposal in Ref. 9|, that we shall refer to as the MK protocol, is restricted to systems 
where the inhomogeneous broadening is provided by the Doppler effect. 

A variant of the MK protocol has been proposed for solids in which the inhomogeneous 
broadening is linked to stochastic variations in atomic transition frequency that depend on 
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an absorbing center's position in a host medium 
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12j. Althou gh there are several proof 
of principle experiments of this absorption-type protocol Ij], 16], all such experiments 
have involved classical input fields. It should be stressed that the large time-bandwidth 
product capacity is lost when the Doppler shift is replaced by a more stochastic source of 
inhomogeneous broadening. In order to keep control of the inhomogeneous phase shift, one 
selects a narrow spectral group of atoms at the beginning. An external field is used to 
spread this initial ensemble over the desired bandwidth. Hence the bandwidth is increased 
at the expense of the available optical density, i.e. at the expense of the capacity to trap 
the optically carried information within the material. 

A key feature in the MK scenario is that the totally absorbed signal is restored without 
amplification. This contrasts with previous photon echo investigations where large retrieval 
efficiency results from strong amplification in an inverted medium [l7, Ij], 19, 20, 21]. There- 
fore, unlike these earlier works, the MK protocol is free from the noise associated with spon- 
taneous and stimulated emission. This is closely related to the fact that the control fields 
do not interact with highly populated states. As a consequence few atoms are promoted to 
the upper electronic level. Another consequence is that, as the control fields interact with 
quasi-empty states, they are neither attenuated nor distorted as they travel through the 
active medium. 

Quite surprisingly, the original MK scheme has not been demonstrated experimentally. 
A possible issue in atomic vapors is the short lifetime of the active optical transition upper 
level. Indeed, to combine a large optical density with the absence of collisions one has to 
work on strong lines with short upper level lifetime. As a consequence, some operations have 
to be carried out on a nanosecond time scale. Specifically, one needs nanosecond 7r-pulses to 
convert optical dipoles into Raman coherences and conversely. In addition, the input signal 
duration is limited to a few nanoseconds. 

In this paper, we propose a Raman variant of the MK scheme that circumvents these 
limitations. Direct excitation of ground state coherence avoids the introduction of rapidly 
decaying quantities, yet retains the other advantages of the MK protocol. Quantum mem- 



ories based on Raman scattering have been proposed in the past |15l . [22j, [23J. However, 
previous proposals did not fully examine the dynamics of the system, focusing on steady- 



state conditions, 



22| or they considered a situation in which the temporal sig nal profile is 



mapped into a spatial distribution of the atomic ground state coherence 15|, |23| . In our 
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scenario, the spectro-temporal features of the signal are stored in the spectral distribution of 
the ground state coherence. The experimental investigation presented in ref. [ijj] is actually 
very close to our situation, but the authors resort to a reversible magnetic field gradient 
to cover the signal spectrum and to reverse the atomic phase, which ultimately leads to a 
spatial mapping of the signal. In our case we use optical pulses to reverse the atomic phase. 

Our objective in this paper is to describe the underlying physics of the Raman protocol, 
a goal that can be achieved within the confines of a theory in which all radiation fields are 
treated classically. The paper is arranged as follows: after presenting a picture of the overall 
process in section II, we develop a theory for each step of the protocol in sections III-V. In 
section VI we discuss the range of applicability of the storage method. 

II. OUTLINE OF THE STORAGE AND RETRIEVAL PROCEDURE 

The protocol consists of the four steps shown schematically in Fig. [U We consider an 
ensemble of four-level atoms. The four levels can be magnetic state sublevels of the same or 
different ground state hyperfine levels. The atoms are prepared initially in state \a) . 

1. Stage 1 

In the first stage, a quasi-monochromatic control field and the input pulse drive Raman 
transitions between levels a and c. Each field is off-resonant for optical excitation of level b, 
but the difference of the field frequencies is close to that of the Raman transition. The con- 
trol and input fields have a relative propagation vector K that leads to a Doppler shift K • v 
associated with the two-photon Raman transition. As a consequence the bandwidth that 
can be absorbed in this Raman process is on the order of 1/Ku, where u is a characteristic 
atomic speed. Of critical importance is that all frequency components of the signal field are 
depleted in an identical fashion (no pulse distortion) if the bandwidth is less than the inho- 
mogeneous width. The medium is optically dense so the signal pulse is totally attenuated. 
In other words, as a result of stimulated Raman scattering, the signal pulse energy is totally 
transferred to the control field. Each atom has a negligibly small population in state |c) and 
the entire population stored in level c is assumed to be small as well, assuming the signal 
pulse is weak. The control field is turned off following the depletion of the signal pulse. In 
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FIG. 1: (Color online) Level scheme and protocol steps, (a) the weak input signal, combined 
with a strong control field, drives the Raman transition a — c. (b) the Doppler phase build-up is 
stopped by conversion of the coherence p ac into p a d- This is accomplished by 7r-pulse excitation of 
the Raman transition c — d. (c) the coherence p ac is recovered with the help of a second Raman 
7r-pulse. (d) a backward propagating control field creates a coherence p a i, which allows for the 
restoration of the signal pulse, propagating in the backward direction. 

contrast to the MK protocol in which an optical coherence is created in the first stage, the 
signal field is transferred directly to a Raman coherence in our protocol that is immune to 
spontaneous emission decay. 

2. Stage 2 

The Raman coherence dephases following excitation as a result of the inhomogeneous 
broadening. As in a photon echo experiment, this dephasing can be reversed by the ap- 
plication of a second pair of pulses. However, if we were to use a two-pulse echo process, 
the Raman coherence excited by the signal would be affected by velocity changing collisions 



over the entire storage time between the two pulses. Instead, if we use a three-pulse echo 
configuration, the effect of collisions during most of the storage time can be suppressed. 
The second pulse pair, in effect, freezes the Doppler phase created by the first pulse pair. 
As in the MK protocol, this phase reversal must be carried out with a Raman pulse that 
leaves the population in level a unchanged. This is a crucial condition to ensure uniform 
illumination by the control fields; if any control field is resonant with a transition originating 
in level a, the field will be strongly absorbed in the optically dense medium and unsuitable 
for this protocol 25[. All these requirements are satisfied if one applies a Raman ir pulse 
between levels c and d, following the excitation of the a — c Raman coherence. The relative 
K vector of the two fields is the same as that in stage 1, so the net effect of the ir pulse is 
to convert the a — c to a — d coherence, while freezing the Doppler phase evolution of this 
Raman coherence, just as in a stimulated photon echo. The net result is that the original 
pulse information is now stored in a Raman coherence that is, for the most part, "protected" 
from the effects of velocity-changing collisions. 



3. Stage 3 

To prepare the system for the retrieval stage, a ir pulse is sent into the medium at some 
later time to restore the a — c coherence. This Raman pulse has its relative K vector reversed 
and prepares the atoms with a spatial Raman coherence that allows for retrieval of the signal 
in stage 4. 



4- Stage 4 

A control pulse is sent into the sample that is identical to the initial control pulse, but 
with its propagation vector reversed. The first three stages produce a Raman coherence that 
allows a field to build up in a direction opposite to that of the input signal pulse. In other 
words, the depletion of the signal field is reversed and the original signal pulse is restored as 
it exits the sample propagating in a direction opposite to that of the original signal pulse. 
In principle, the pulse can be restored with close to 100% fidelity. In the context of creating 
a functional quantum memory device, splitting the phase reversal in two steps (stages 2 and 
3) reduces the waiting time between the read-out pulse and the emission of the restored 



6 



signal, making the retrieved data available more rapidly after the read-out decision. 
We now describe in detail how each of these steps can be achieved. 

III. STORAGE STEP 

Storage corresponds to the mapping of the optically carried information into atomic 
Raman coherence, accompanied by attenuation of the input signal field. 

A. Buildup of a Raman atomic superposition state 

The input signal is depleted by stimulated Raman scattering in a three-level A-system. 
The weak pulse to be stored, combined with the control field, resonantly excites the Raman 
transition \a) — |c). Both fields are tuned off resonance from the optical transition to upper 
level | b). Storage is performed into the superposition of the ground substates \a) and |c). 
Since the atoms are prepared initially in state \a), the medium is transparent to the control 
field that uniformly illuminates all the active atoms. Assuming that the control field is 
constant during the signal pulse, we write the electric field of the control field as the plane 
wave 

E 2 (r,t)=A 2 e^- r - i ^ t + c.c (1) 

where the envelope A 2 is a time- and space-independent parameter. The control field wave 
vector and frequency are denoted k2 and u>2, respectively. 

The signal pulse Rayleigh range is assumed to be much larger than the storage material 
length L to insure that its diameter does not vary significantly as it propagates in the 
medium. The electric field of the signal field can be expressed as 

£i(r, t) = Ai(r, t)e ikl - r ~ luJlt + c.c, (2) 

where A\(r,t) is the envelope, ki the propagation vector, and uj\ the carrier frequency of 
this field. The spatial dependence of A\(r,t) reflects the radial distribution of the field and 
its attenuation along direction ki. When L/c is not much smaller than the pulse duration, 
retardation also contributes to the field envelope spatial dependence. The Rabi frequencies 
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associated with the signal and control fields are denoted by 

Q 1 (r,t) = -ii ha A 1 (r,t)/h; (3a) 
^2 = -PbcA 2 /h, (3b) 

where pba and pbc are optical dipole moment matrix elements. It is assumed that k\ ~ k 2 = k. 

In perturbation theory with the population of level a set equal to unity, the coupled 
equations for the optical and Raman coherences are: 

p ah , m = (iAi - lab )p ah . m + zfij [r m (t),t] e" ikl - r -W + zp ac;m n; e - lk2 - rm(t) (4a) 
p ac , m = [i(Ai - A 2 ) - 7 ac ]p ac;m + zpafe;mn 2 e ik2 - r " l(i) (4b) 

where Ai and A 2 are the atom-field detunings for each optical transition, 7^ is the decay 
rate for the a — b coherence, 7 ac is the decay rate for the a — c coherence, p a b;m = Pab;m^~ tuJlt 
and p aC ;m = Pa^m^^ 1 "^ ■ These equations give the time evolution of the density matrix 
elements for atom m, located at r m (t) at time t. The manner in which the spatial phases of 
the field are imprinted on the atoms is readily apparent in Eqs. (pD). Under the assumption 
that Qi, Qi dQi/dt, 7^, ku « Ax, where u is the most probable atomic speed, the optical 
coherence adiabatically follows the field variations and can be written as: 

Pab-, m = -fij [r m (t), t] e'^^ / A, - p ac . m Q* 2 e- ik2 - r ^ / A 1 . (5) 
Substituting this expression into Eq. |4bl one obtains 



in 2 | 2 



Ax - A 2 - i^J- I 



, nj [r m (t),t] n 2 e tK - r -w 

Pac;m ~ % T (6) 

^1 



where 

K = k 2 -k!. 

The detuning Ai — A 2 can be adjusted to cancel the light shift |n 2 | 2 j A\. Finally, the Raman 
coherence can be expressed as: 

Pac;m(t) = f ^ [v m (t'),t'} e**^ (7) 

where it has been assumed that any decay of p ac -m during the signal pulse can be neglected. 
The density matrix element can be expressed in terms of the atom position at time t. Indeed, 
if collisions do not change the atomic velocity v m during the signal pulse, the position at t' 
can be expressed as r m (t') = r m (t) — v m (t — t'), so that: 

Pac-At) = -i^e*-*"® f dt'Q{ [v m (t) - v m (t - t'),t'} e^^') (8) 

^1 J-oo 
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B. Stimulated Raman scattering of the input signal 



It is assumed that, in the absence of the control field, the scattering of the probe field 
is negligible owing to the large detuning Ai. As a consequence, the contribution to the 
polarization resulting from the first term in Eq. (j5J) can be neglected. On the other hand, 
the control field intensity is sufficiently large to allow for significant stimulated Raman 
scattering, resulting in a loss of signal field intensity as the signal field propagates in the 
medium. The Raman contribution to p a b;m is given by the second term on the right hand 
side of Eq. (jSj). To determine the modification of the signal field, we need to calculate the 
polarization associated with the p a b- m coherence. We write the macroscopic polarization as 

P(r, t) = P+(r, t) e ik i- r -^* + P*( r , t) e -(.*i-'-«>**) 

In going over to a macroscopic polarization, we assume that, at any position r in the 
medium, one can define a slice of thickness I << 2ir/k in the ki direction, containing many 
atoms. The macroscopic polarization at (r, t) is obtained by combining the contributions 
from all the atoms within the slice. Those atoms satisfy the condition ki.s m (t) < 1/2, 
where ki is a unit vector in the ki direction and s m (t) = r — r m (t). Therefore the positive 
frequency component P + (r,t) is given by 

P + (r,t) = ^e-^- r £ p ba -, m (r,t), (9) 

m 

|ki.s m (f)|<Z/2 

where SV represents the slice volume. Combining Eqs. (151), (Q, and ([5a]) . we find 

P +M =i WMZ ^ / dt'Mr-Vmit-t'U}^^. (10) 

|ki.s m (i)|<//2 

As was noted above, the contribution to p& a;m (r, t) from the first term in Eq. (jSJ) has been 
neglected since it adiabatically follows the field and vanishes for times greater than the pulse 
duration. 

The atoms are uniformly distributed in space, with density N, and their normalized 
velocity distribution is represented by W(v). Replacing the discrete sum by an integral, 
according to 

^ £ - N I d 3 vW(v), 

m J 
|k 1 .s m (t)|<Z/2 
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enables us to transform Eq. (TTOT) into 

I 1 2 I 1 2 t 

P + (r ; t) = i ^ ab [^ N [ d 3 vW(v) [ dt'Ax [r - v(t - t'),t'] ^ {t ^\ (11) 

f^A 1 J J-oo 

Provided the input signal spectral width 5 S is smaller than Ku, A\ [r — v(t — t'),t'] can be 
taken out of the integral over t' and evaluated at t' = t. In this limit the polarization P+(r, t) 
reduces to: 

I . |2 |Q |2 p poo 

P + (r,t) = i ™ ' 2| At(r,t)N / dvW{v) / rfre^ r (12a) 

J Jo 



i N7r\fi ab \ 2 W{0)k\Q 2 \ 2 



At{r,t) (12b) 



k H KAj 

where W(v) represents the one-dimensional velocity distribution. 

This is the key result of this section. Owing to the large inhomogeneous width, the 
polarization is proportional to the field amplitude and depends locally on this amplitude. In 
other words, the polarization does not depend on the value of the field amplitude at earlier 
times as it would in the case of homogeneous broadening. The electric susceptibility xr 
is defined by P+(r, t) = eoXif,^i( r ^) with the intensity absorption coefficient o>r given by 
a/; = klm(xR), which leads to 



Q R = r/A2 Q ( 13 ) 



|2 

KAf 

where ao represents the linear absorption coefficient on the inhomogeneously broadened 
transition a — b. Provided 5 S < Ku, the signal propagates without distortion, which implies 
that all the signal spectral components are uniformly attenuated and stored in the atomic 
ensemble. The spectral components are mapped into Raman coherence in atoms whose 
velocities span an interval of order 5 S /K along the direction K = k 2 — kj. The parameter 
K should be adjusted in such a way that S s /Ku is larger than unity to provide a sufficiently 
large bandwidth to store the signal pulse, but not too large, since the signal depletion and 
Raman storage varies inversely with K. 

IV. FREEZING THE DOPPLER PHASE 

The input signal illuminates the storage medium during a time interval centered at t\. 
The control field is turned off following the signal pulse. At some later time, the Raman 
coherence, given by Eq. ([8]), can be expressed as: 

Pac;m(t) = -z e iK ^ t ^ iK - V - + ^ t -^n m (14) 
10 
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FIG. 2: (Color online) Schematic representation of the timing sequence for the entire protocol, 
showing the relative duration of the different pulses, and the various atomic quantities involved at 
the different stages. 



where 

q roo 

K m = -^ drOi [r m (ti),h + T]e tK - v - T (15) 

and it was assumed that [r — VT p ,t] ~ Q\ [r,t], where r p is the signal pulse duration. 
There is a build-up of Doppler phase associated with the Raman coherence that grows 
linearly as a function of time following the interaction with the signal pulse. Although 
this phase could be reversed at a later time, it is best to nip it in the bud to prevent any 
deterioration from velocity-changing collisions. To accomplish this task one can send in a 
Raman 7r pulse having the same K vector that transfers the amplitude from state |c) to an 
auxiliary level \d), or, equivalently, converts coherence p ac into p ad . 

Let Q 3 (t) and f2 4 (t) denote the Rabi frequencies on transitions |c) — \b) and \b) — \d), re- 
spectively. The radial extension of these fields is assumed to be larger than that of the signal 
pulse. Moreover, these fields are not attenuated as they propagate through the storage ma- 
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terial since interact with quasi-empty levels. Hence their Rabi frequency space dependence 
can be omitted. Both fields are detuned by the same amount A from resonance with their 
respective assigned transition. Therefore the |c) — \d) two-photon transition is resonantly 
excited. The Rabi frequency of the equivalent two-level system is given by Q 3 (t)Ql(t)/A. 



The two-level approximation holds provided 



f2 3i 4(A ± 6 S ) << 1, where Q(A) represents 



the time-to-frequency Fourier transform of Q(t). Of course one must take care that none 
of those fields can excite state \a) , where all the atomic population is concentrated. With 
these assumptions, one finds that p ac and p ad evolve according to 

Pad;m = e Pac;m ~ ladPad;m (16a) 

2 

6 ^ ^ ^ Pad;m ''lacPac-^mi (16b) 



where p ad . m = p ad ; m e Ui)t , <p m {t) = (k 3 - k 4 ).r m (t) and, for simplicity, we have set 
f2 3 (t) = (t) = Q n (t) and neglected a light shift that is the same for levels c and d. The 
pulses, having duration T n , are applied at a time centered about t = t 2 and their duration 
is assumed to be sufficiently short, SgT^ << 1, to resonantly excite all the atoms that were 
excited by the signal pulse (recall that atoms in the velocity range |K ■ v| < \5 S \ are excited 
in stage 1). The time-diagram of the whole protocol is displayed in Fig. [2J showing the 
relative duration of the different pulses. 

In terms of the coherence components at time t 2 just before the Raman ir pulse is applied, 
Eq. (TTBT) can be solved, for t > t 2 , as: 

~Pac-, m (t) = e-^-^ cos f Ij PacA^) + le 1 ^-^-^ sin (®j p« d;m (*2~) (17a) 
Pad;m(t) = ze^^ 2 )"^" 42 ) sin p ac;m (t 2 ") + e^" 2 ) cos p ad;m (t 2 -), (17b) 

where 6 = 2 / dt |Or(t)| /A, assuming that 7 ac 7V , Tad^ « 1- Total conversion from p ac to 
p a d requires that = n. With initial conditions given by Eq. ffT4l) and p a d;m(t 2 ) = one 
finds that 

padm(t) — e^ k4_k3+K ^' rm( '* 2 - )_7ad ^ _ * 2 ' ) e _ ^ K ' Vm+7ac ' )tl2 7?. m (18) 

where is the time interval between the j'th and ith pairs of pulses. With k 3 — = K, 
there is no build-up of Doppler phase following this second Raman pulse 26|] - the Doppler 
phase has been frozen. 
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V. SIGNAL RECOVERY 



To prepare for signal retrieval at some later time, one sends in another Raman it pulse, 
centered at time t — t^, to restore the p aC ;m coherence. This Raman pulse consists of two 
fields having propagation vectors and kg that drive the \d) — \b) and \b) — \c) transitions, 
respectively. The atom-field dynamics is again described by Eqs. fll6l) . with <f) m (t) replaced 
by (j)' m (t) = (k' 3 — \t' A ).r m (t). To reverse the Doppler phase acquired in the time interval 
t\2, we choose the propagation vectors such that k' 3 — k^ = — K, which leads to a Raman 
coherence for t > t 3 given by 

X — ? V0Jn(*3)-2K.V m tl2-7a C (i-i3+il2)-7ad*23<7? 

_ ^ e -*K.r m (t)+iK.v m (t-t 3 )-iK.v m ti2-7ac(i-i3+ii2)-7ad*23-^ 

To restore the signal one directs a quasi-monochromatic control field having frequency u 2 
into the medium in the backward direction, with k 2 ' = — k2. This field should be switched on 
somewhat before time t = t 3 + t u , with the same Rabi frequency fl 2 as the initial coupling 
beam. The restored field can be written as 

E{r,t) =^(r,t)e ik " r - iWl * + c.c. 

We now argue that, only if k" = — kx, can a phase matched signal be generated. To see 
this, we return to Eq. (fT9l) and analyze the phase factor in this equation. Although lZ m 
contains a Doppler phase factor [see Eq. (jl5p ]. it is of order Kut, where r is the signal 
pulse duration, and is small compared to the other Doppler phases appearing in Eq. (FTPl) . 
since we assume that t\2 3> r. When the control field is sent into the medium, the Raman 
coherence ( TT9]) gives rise to an optical coherence on the a — b transition that is responsible 
for the generation of the restored field. ^From Eq. (jSJ), one can deduce that the optical 
coherence p a b;m varies as 

e ik 2 .r m (t)-iK.r m (t)+iK.v m (t-t 3 )-iK.v m ti2 
_ e iki.r m (t)+iK.v m (t— *3— tia) 

This expression implies that a phase matched signal can propagate only in the k" = — ki 
direction and that this signal can be nonvanishing (owing to the average over atomic veloc- 
ities) only for times t pa t 3 + 1 12 . Thus, in what follows we assume that the restored signal 
field has propagation vector k" = — k^ 

13 



The existence of the restored field must be taken into account in the expression for the 
Raman coherence. Indeed, this field also combines with the control field to drive Raman 
transitions between levels a and c. Adding this contribution to Eq. (FTU]) . one obtains the 
Raman coherence 

X (+\ — ,;„-iK-.r m (t)+iK.v m (t-t3-t 1 2)-2-y ac t 1 2--y a dt23'Tp 

.^ab^2 f l 1,1 , r _ u ,s ,n AK".r m (t') ^ 



Substituting this equation into Eq. (jSj) one arrives at an expression for the optical coher- 
ence, 

|2 



/^ qfe|^2| a ,ki, fa (t) 



X | e -7a C (t-*3-*12)-7a d *23 f dr ^* ( f 1 ) ? ^ + r ] e «K.V m (t-t 3 -t 12 +r) (g^ 

>/ —00 / 
that can be used to calculate the macroscopic polarization density, which is written as 

P s (r,t) = P +s {r,t)e~ l]il - r - luJlt + P* a (r, t) e -(- ik i- r -^*) j 

incorporating the fact that the signal is phase- matched in the — ki direction only. The 
polarization density is comprised of two components. One of them represents the source 
term that gives rise to the restored signal. This contribution originates from the first term 
on the right hand side of Eq. (l2Ti) and is given by 

|2 



pWfr +\ - - l/^frl l^2| N_ - 7ae ft-t 3 +t 12 )-7 arf t 2 3 

r +a [r,t)- i e 

/POO 
dvW(v) / drA! [r(ti), h + r] e -iK.v(t-t 3 -ti2+r) 
•/ —00 



where A\ [r(ti),t'] ~ Ai[r — vti 3 ,t']. The polarization density depends on the field expe- 
rienced by the participating atoms at their positions when they interacted with the signal 
field. 

For the field to be restored, the condition A\ [r m (ti), t'] ~ A\ [r, t'] must hold. This 
condition is valid provided the distance travelled by the atoms in time is much smaller 
than the beam diameter, the spatial width of the input signal envelope, and the absorption 
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length a R . Since 5 S « Ku, one can take Ai [r, t'\ out of the integral over t' and evaluate 
it at time t' = t\ — (t — — £12), which leads to 



(23) 



/poo 
dvW{v) / dTe~ iKvT 
J — oo 

Using Eqs. (fT2j) and ({TBI) , one can re-express this in terms of «s as 

P«(r, f) = -2*6 ^A [r, ^-(t-ts- fia)] e^*--^ (24) 

The amplitude Ai [r, ti — (t — £3 — £12)] is nonvanishing only for t ~ £3 + £12, reflecting the 
fact that Doppler rephasing occurs only for such times. We have used this fact to set 

g — ~fac(t — *3— tl2) po g — 27ac*12 

The input signal was attenuated as it traveled through the storage medium. Taking into 
account the propagation and attenuation of the input field A\ [r, t] , we can write Pj. (r, t) 
in terms of the input field at z = as 

P«(z, t ) = -2ie ^Ai [0, h-(t-t 3 - t 12 ) - z/c] e -^z/\-2 lac t, 2 - lad t 2 ^ (25) 

the z axis being directed along ki. 

The contribution to the polarization density from the second term on the right hand 
side of Eq. (1211) corresponds simply to depletion of the restored field as a result of Raman 
transitions from level a to c . This component can be written as 

P®(z t t)=ie ^A(z,t). (26) 

Finally, the restored field is a solution of the linearized Maxwell equation: 

- d,A(z,t) + -d t A(z,t) = <A \p$(z,t) + Pl 2 ](z,t)} (27) 

Making the change of variables z' — z, t' — t + z/c, and combining Eqs. (125j) and ( 1271) . one 
obtains 

d z ,A'(z', t') = -a R A x [0, h - (f - t 3 - t 12 )\ e -^/ 2 e- 2 ^-^ + ^a^V, 0, (28) 

where A'(z', t') = A(z', t' — z 1 /c). Solving this equation with the initial condition *4/(L, t') = 
0, we find 

A(z, t) = Ax [z, h-(t-t 3 - t 12 )} e- a T«*"-T-«'*» (1 - e - Q *( L -^) . (29) 
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For olrL 3> 1, the pulse is totally restored at time t = t s + tu , neglecting decay of the 
Raman coherences. The latter equation represents the main result of the paper, showing 
the absence of distortion of the retrieved signal and its variation as a function of the optical 
depth ocrL. 

One can also verify that the atoms return to their initial state as the signal field is 
restored. Substituting Eq. (129]) into Eq. (|20|) and using Eq. (Tl5l) . one obtains 

% (+\ _ - /fnjA? -iK.r m (t)+iK.v m (t-t 3 -ti2)-7ac(2t23+ti2) 



X 



( r dt'AKz, t 'y K - v ^-^ - (i - e - aR{L ~ z) ) r dt'AKz, t 'y K - v m(f-ti) \ 

\J-00 J tl-(t-t 3 -tl2) J 



(30) 

For times (t — t 3 — t 12 ) > the lower limit on the second integral can be replaced by — oo 
and this term cancels the first integral. In other words, the coherence vanishes once when 
the signal field is restored, provided a^L 3> 1 . 

Following the emission of the restored signal pulse, there remains in the medium a Raman 
coherence given by 

p ac .(t) = ^ ab ^ 2 e -iK.r ro (t)+tK.v m (t-t 3 -ti2)-7 ae (t33+2ti2) e ~a R {L-z) f dt'A*(z f'V^'^C*'-*!) 

(31) 

Comparing with expression with Eq. (151) , and noting that the population in state \a) is 
approximately equal to unity, one can conclude that the population left in state |c) is 



PccAih +t 12 ) + ] = \pac;m(tt)\ [l - e~^ tM ^] + \p ac , m [(t 3 + t 12 ) + ] [ (32a) 

= \p aC ; m (tt)\ 2 [1 - e- 2 ^+ 2 *-) (1 - e-***^*))] . (32b) 



Finally, summing over z, one finds the fraction rj of the initially excited atoms that is left 
in state |c) is given by 

r) = 1 - e - 2 T-( t23+2 * 12 ) (1 - e~ aRL ) . (33) 

If decay of the Raman coherence decay can be neglected, rj = e~ aaL . During the storage 
process, a fraction t] in = (1 — e~ aRL ) of the incoming signal radiation is used to excite the 
atoms to level c , the remaining signal passes through without being scattered. From the 
part that is stored, a fraction is lost at retrieval, even in the absence of relaxation. Indeed 
the restored field is (1 — e~ a ° L ) times smaller than the incoming one, according to Eq. 



1(3 



Therefore one recovers a fraction W ou t/Wi n = (1 — e _QflL ) 2 of the incoming energy. The 
difference r\i n — W ou t/Wi n = e~ aRL corresponds to the fraction rj of the excited atoms that 
is left in level c. To summarize, with a finite length material, information is lost in equal 
amounts at storage and retrieval, the storage loss associated with an incomplete depletion 
of the input field and the retrieval loss associated with excited state population p cc that is 
left in the medium. 

VI. DISCUSSION 

We have shown that nearly 100% recovery efficiency can be reached with the proposed 
quantum storage Raman scheme. The memory bandwidth, determined by the Raman tran- 
sition Doppler broadening, can be tuned continuously from to 2ku. Since the atoms are 
not excited to the upper electronic level, the control pulse duration is not limited by the 
upper level lifetime, but only by the inverse signal bandwidth. 

Our motivation in this work was to circumvent the time scale conditions imposed by the 
short upper level lifetime, since it is not easy to produce large area, large waist coherent 
pulses on the nanosecond time scale. The pulses to be considered here, whether input signal 
or control pulses, will most likely last for tens or hundreds of nanoseconds. This corresponds 
to a spectral width much smaller than the Doppler width ku, typically of order lO 9 ^" 1 or 
higher. The available bandwidth of the Raman scheme is determined by Ku, where K/k 
can be expressed in terms of the angle 9 = (k l7 k 2 ) as K/k = 2#sin(6 l /2). According to 
Eq. (129]) . the Raman optical depth a^L should be as large as possible. Since a^L is 
inversely proportional to K/k, the value of K/k should be matched to the bandwidth of 
the input pulse. For an input field spectral width of order 10 MHz, we are led to a nearly 
co-propagating configuration with 9 of order 10 mrad. 

The memory lifetime is limited by the atomic motion. As noticed in Section IV, the 
distance travelled by the atoms during the entire process should be much smaller than the 
beam diameter and the absorption length a^ 1 . With a typical average speed of a few lOOm/s, 
a memory lifetime of a few tens of microseconds limits the sample size to a few centimeters. 

A critical issue in a multilevel system is the ability to selectively drive the target tran- 
sitions. Field polarization can be used to provide the selectivity if the angle between the 
fields is small (or close to 7r); in this limit, one can use circularly polarized fields, as shown 
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FIG. 3: (Color online) Example of a possible excitation scheme. The Raman transitions connect 
the Zeeman sublevels of F = 1, F = 2 hyperfine levels. Each Raman transition involves a <r+ and 
a <t_ polarized beam. 

in Fig. [3j The ground state manifold consists of F — 1 and F = 2 hyperfine states and the 
A-systems of the Raman protocol involve magnetic substates of these levels. Since the angle 
between the beams is assumed to be on the order 10 mrad or so, the fields can taken to be 
cross-polarized, in first approximation. 

Optical pumping by the strong control field determines the uncoupled initial state (F = 1, 
m — — 1). In step (a), fields 1 and 2 create a Raman coherence where the notation 

is PF,m:F>,m>- The Doppler dephasing is then stopped by a 7r-pulse, composed of two beams 
propagating along k 3 and k 4 in such a way that k 4 — k 3 = ki — k 2 [step (b)]. These fields 
convert the coherence pi-\a t \ into pi,_i : 2,-i- The process in which fields 4 and 3 drive 
Raman transitions between states |1,— 1) and |1,1) is suppressed, provided the hyperfine 
frequency separation is much larger the inverse pulse duration. In step (c) a Raman n pulse 
having k 4 — k 3 = — (k 4 — k 3 ) restores the /Oi_i : i,i coherence and prepares the system in such 
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a fashion that the subsequent application of a control pulse having propagation vector — k2 
restores the input signal field propagating in the — ki direction at time t — ts + tn- 

VII. CONCLUSION 

We have proposed a Raman based quantum memory scenario that operates with close to 
100% efficiency. This protocol fits well to intermediate time scales, with signal duration of 
order 100ns, a time range that is well adapted to experiments based on high spectral purity 
continuous wave laser sources. We have shown that the spectral components are stored 
in different atomic velocity classes and have explained how to optimize the optical depth 
depending of the needed storage bandwidth. Since we have assumed that the signal field is 
weak and since spontaneous emission is negligible in this protocol, we expect the results to 
be unchanged if the classical input field is replaced by a quantized, pulsed radiation field. 
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